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Abstract 

The scaling limit of the two-dimensional Ising model in the plane of temperature and mag- 
netic field defines a field theory which provides the simplest illustration of non-trivial phe- 
nomena such as spontaneous symmetry breaking and confinement. Here we discuss how Ising 
field theory also gives the simplest model for particle decay. The decay widths computed 
in this theory provide the obvious test ground for the numerical methods designed to study 
unstable particles in quantum field theories discretized on a lattice. 



1 To appear in the proceedings of the XVth International Congress on Mathematical Physics, Rio de Janeiro, 
6-11 August 2006. 



1 Ising field theory 



Quantum field theory provides the natural tool for the characterization of universality classes of 
critical behavior in statistical mechanics. While the general ideas based on the renormalization 
group apply to any dimension (see e.g. [JJ), the two-dimensional case acquired in the last decades 
a very special status. Indeed, after the exact description of critical points was made possible by 
the solution of conformal field theories [2| , it appeared that also specific directions in the scaling 
region of two-dimensional statistical systems can be described exactly [3]. Additional insight 
then comes from perturbation theory around these integrable directions [JJ. 

The two-dimensional Ising model plays a basic role in the theory of critical phenomena since 
when Onsager computed its free energy and provided the first exact description of a second order 
phase transition [5] . Its scaling limit in the plane of temperature and magnetic field defines a field 
theory - the Ising field theory - which provides the simplest example of non-trivial phenomena 
such as spontaneous symmetry breaking and confinement [6]. Here we will discuss how Ising 
field theory also yields the simplest model for particle decay [7]. 

The Ising model is defined on a lattice by the reduced Hamiltonian 

E = -— ^(TiCTj ~ H ^CTi , <Tj = ±l (I) 



T 



so that the partition function is Z = Y2{cn} e ■ On a regular lattice in more than one dimension, 
the model undergoes, for a critical value T c of the temperature and for vanishing magnetic field 
H, a second order phase transition associated to the spontaneous breakdown of spin reversal 
symmetry. 

In two dimensions the scaling limit of ([1]) is described by the Ising field theory with action 

A = Ac ft ~ r J d 2 xe(x) - h j d 2 xa{x) . (2) 

Here Ac ft is the action of the simplest reflection-positive conformal field theory in two dimen- 
sions, which corresponds to the Ising critical point [2]. The spin operator cr(x) with scaling 
dimension X a = 1/8 and the energy operator e(x) with scaling dimension X e = 1 are, together 
with the identity, the only relevant operators present in this conformal theory. The couplings h 
and r account for the magnetic field and the deviation from critical temperature, respectively. 

The field theory ([2]) describes a family of renormalization group trajectories flowing out of 
the critical point located at h = r = (Fig. 1). Since the coupling conjugated to an operator <1> 
has the dimension of a mass to the power 2 — X$, the combination 

is dimensionless and can be used to label the trajectories. In particular, the low- and high- 
temperature phases at h = and the critical isotherm r = correspond to rj = — oo, +oo, 0, 
respectively. 
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Figure 1: The Ising field theory ([2|) describes a one-parameter family of renormalization group 
trajectories (labelled by rj) flowing out of the critical point located at r = h = 0. 

2 Evolution of the mass spectrum 

Under analytic continuation to imaginary time the euclidean field theory ([2]) defines a 1+1- 
dimensional relativistic theory allowing for a particle interpretation. 

It is well known that the action ([2]) with h = describes a free neutral fermion with mass 
m ~ |t|. While in the disordered phase r > this fermionic particle corresponds to ordinary 
spin excitations, in the spontaneously broken phase r < it describes the kinks interpolating 
between the two degenerate ground states of the system (Fig. 2a). In the euclidean interpretation 
the space-time trajectories of the kinks correspond to the domain walls separating regions with 
opposite magnetization. 

A small magnetic field switched on at r < breaks explicitely the spin reversal symmetry 
and removes the degeneracy of the two ground states (Fig. 2b). To first order in h the energy 
density difference between the two vacua is 

A£ ~ 2h {a) , (4) 

where (a) is the spontaneous magnetization at r = 0. With the symmetry broken the kinks are 
no longer stable excitations. An antikink-kink pair, which was a two-particle asymptotic state of 
the theory at h = 0, now encloses a region where the system sits on the false vacuum (Fig. 2b). 
The need to minimize this region induces an attractive potential 

V(R)~A£R (5) 

(R is the distance between the walls) which confines the kinks and leaves in the spectrum of the 
theory only a string of antikink-kink bound states A n , n = 1, 2, . . ., whose masses 

m n = 2m+ {A£) ^ Zn , h^O (6) 

are obtained from the Schrodinger equation with the potential ©. The z n in ([6]) are positive 
numbers determined by the zeros of the Airy function, Ai(— z n ) = 0. This non-relativistic 
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Figure 2: The free energy and the kink excitations in the spontaneously broken phase (a). A 
small magnetic field removes the degeneracy of the ground states and confines the kinks (b). 

approximation is exact in the limit h — > in which m n — 2m — > 0. The spectrum ([6]) was first 
obtained in [6] from the study of the analytic structure in momentum space of the spin-spin 
correlation function for small magnetic field. Relativistic corrections to © have been obtained 
more recentlj^] in [5] . 

The particles A n with mass larger than twice the lightest mass mi are unstable. It was 
conjectured in [6] that the number of stable particles descreases as rj increases from — oo, until 
only Ai is left in the spectrum of asymptotic particles as r\ — > +oo. The particle A\ would 
then be the free fermion of the theory at r\ = +oo. According to this scenario, for any n > 1 
there should exist a value r] n for which m n crosses the decay threshold 2mi, so that the particle 
A n becomes unstable for rj > rj n . The natural expectation is that the values % decrease as n 
increases, in such a way that the r-h plane is divided into sectors with a different number of 
stable particles as qualitatively shown in Fig. 3. The trajectories corresponding to the values r\ n 
are expected to densely fill the plane in the limit i] — ► — oo. 

This pattern has been confirmed by numerical investigations of the spectrum of the field 
theory ([2]) for all values of rj [H [HI [9]. For r\ — > — oo the particles A n with n large can be studied 
within the semiclassical approximation and their decay widths have been obtained in |1U[ |2] ■ 

3 Particle decay off the critical isotherm 

The critical isotherm r/ = r = must lie within the sector in which the theory has, generically, 
three stable particles (Fig. 3). To understand this we must recall that A. Zamolodchikov showed 
in [3] that the theory (|2|) with r = is integrable and computed its exact S-matrix. He found 
that the spectrum along this integrable trajectory consists of eight stable particles A\,... ,A% 

1 See also [9] which appeared after this talk was given. 
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Figure 3: Expected evolution of the mass spectrum as a function of r\. In the sector in between 
r\ n and ?7 n +i the theory possesses n stable particles (numbers in the circles). 

with masses 
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(7) 



A peculiarity of this spectrum is that only the lightest three particles lie below the lowest 
decay threshold 2mi. The remaining five have the phase space to decay and certainly are not 
prevented to do so by internal symmetries (the magnetic field leaves no internal symmetry in 
the Ising model). It is easy to see that, while there is nothing wrong with the stability of the 
particles above threshold along this integrable trajectory, they must necessarily decay as soon 
as a deviation, however small, from the critical temperature breaks integrability [7j. Figure 4 
shows the bound state poles and the unitarity cuts of the elastic scattering amplitudes Su and 
512 m the complex plane of the relativistic invariant s (square of the center of mass energy) . We 
know from [3] that at r = the scattering channel A%A\ produces the first three particles as 
bound states (Fig. 4a), while the channel A1A2 produces the first four (Fig. 4b). The absence 
of inelastic scattering in integrable theories allows only for the unitarity cut associated to the 
elastic processes. When integrability is broken (i.e. as soon as we move away from r = 0), 
however, the inelastic channels and the associated unitarity cuts open up. In particular, the 
process ^1^2 — ► ^4i^4i acquires a non-zero amplitude, so that the threshold located at s = Am\ 
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becomes the lowest one also in the A1A2 scattering channel (Fig. 4c). Since the pole associated to 
is located above this threshold, it can no longer remain on the real axis, which in that region 
is now occupied by the new cut. The position of the pole must then develop an imaginary part 
which, according to the general requirements for unstable particles [11], is negative and brings 
the pole through the cut onto the unphysical region of the Riemann surface. The other particles 
above threshold, which appear as bound states in other amplitudes at r = 0, decay through a 
similar mechanism. We see then that, beacuse of integrability, the trajectory 77 = corresponds 
to an isolated case with eight stable particles inside a range of values of 77 in which only the 
particles A\, A2, A3 are stable. 

These decay processes associated to integrability breaking can be studied analytically through 
the form factor perturbation theory around integrable models [3]. Indeed, if the action of the 
perturbed integrable theory is 



Antegrable + A / d 2 x*(x), 



(8) 



the perturbative series in A can be expressed in terms of the matrix elements of the perturbing 
operator ^ on the asymptotic particle states (Fig. 5). These matrix elements can be com- 
puted exactly in the unperturbed, integrable theory exploiting analyticity constraints [T21 CE2] 
supplemented by operator-dependent asymptotic conditions at high energies [T4" l [15| [16] . 

For our present purposes we then look at the action ([2|) as the integrable trajectory r = 
perturbed by the energy operator e(x) and must compute corrections in r. The matrix elements 
of an operator <fr(x) in the unperturbed theory can all be related to the form factors^] 



ai ...an 



(o\m\^(0i)...A an (e n )), 



(9) 



where |0) is the vacuum state and the asymptotic states are built in terms of the eight particles 
which are stable at r = 0. 

To lowest order in r the corrections to the real and imaginary parts of the masses come from 
the perturbative terms in Fig. 6 and are given by [I] 



Im mi 



- 2 £ 

a<b , m a +rrn,<.m 



1S„ 



-2r/ c , 

I fcab I 



m c m a 



sinh 9a 



4,5 



with 



Qic ab ) - g dgterxningd by energy-momentum conservation at the vertices in Fig. 6b. 



f c = F^(ilT,0), 

f cab = F! ab ^,9^ b \6l cab) ) 



(10) 
(11) 

(12) 
(13) 



2 The energy and momentum of the particles are parameterized in terms of rapidities as {p°,p 1 ) = 
(m a cosh 8, m a sinh 6) . 
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Figure 4: Poles and unitarity cuts for the elastic scattering amplitudes Su and S\2 in the 
integrable case r = 0, (a) and (b), respectively, and for r slightly different from zero (c). In 
(c) the particle A4 became unstable and the associated pole moved through the cut into the 
unphysical region. 
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Figure 5: Perturbative expansion for a scattering amplitude in the theory ([8]). The matrix 
elements of the perturbing operator ^ can be computed exactly in the unperturbed, integrable 
theory. 
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Figure 6: Terms determining the leading corrections to the real (a) and imaginary (b) parts of 
the masses in Ising field theory at small r. For c > 5 also diagrams with more than two particles 
in the intermediate state contribute to the imaginary part. 





a 


e 


Fi 


-0.640902.. 


-3.706584.. 


F 2 


0.338674.. 


3.422288.. 


F 3 


-0.186628.. 


-2.384334.. 


Fa 


0.142771.. 


2.268406.. 


F 5 


0.060326.. 


1.213383.. 


F 6 


-0.043389.. 


-0.961764.. 


F 7 


0.016425.. 


0.452303.. 


F s 


-0.003036.. 


-0.105848.. 



Table 1: One-particle form factors for the operators a and e in Ising field theory at r = 
|141 117] . The rescaling implied by the notation <!> = ensures that the results in the table 

are universal. 

The available decay channels for the particles above threshold are determined by the spec- 
trum (|7|): 



A 4 - 


+ A X A X 






■+ AiAx, 


AiA 2 


A G - 


■+ A X A X , 


A X A 2 , A X A 3 , A 1 A l A l 



and similarly for Aj and A$. For c > 5 the sum in (jlip must be completed including the 
contributions of the decay channels with more than two particles in the final state. 

One- and two-particle form factors for Ising field theory at r = have been computed in 
[14:\ [TT] (see [18] for a review) . Table 1 contains the complete list of one-particle matrix elements 
for the relevant operators. The results for the lightest particle are compared in Table 2 with 
the result of numerical diagonalization of the transfer matrix on the lattice [19] . Three-particle 
form factors have been computed in [7] in order to determine the imaginary parts f j 1 1 j> . 
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Field theory 


Lattice 


Ff 

n 


-0.640902.. 
-3.70658.. 


-0.6408(3) 
-3.707(7) 



Table 2: Lightest-particle form factors at r = from integrable quantum field theory [141 \T7\ 
and from numerical diagonalization of the transfer matrix |19j . 



Figure 7: Evolution of the mass spectrum as predicted by (I16p as r increases from 0. Ai and A% 
approach the threshold and will decay for positive values of r; A4 has already become unstable 
for r < and moves further away from the threshold. 



The available results relevant for the leading mass corrections are 



h = (-17.8933..)(e) 
h = (-24.9467..) (e) 
h = (-53.6799..) (e) 
U = (-49.3206..) (e) 



I/411I 
I/511 1 
I/512 1 



(36.73044..) |(e) 
(19.16275..) |(e) 
(11.2183..) I (e}| 



where (e) is taken at r = 0. The ratios 



5 Re m 2 n 
hm — — 

??->o 55 



/o 



(14) 



where <5£ is the variation of the vacuum energy density, are completely universal and particularly 
easy to check numerically. These predictions for the variations of the real part of the masses of 
the first four particles have been confirmed numerically both in the continuum [U [9] and on the 
lattice [20]. Using the result [21] 



(15) 



(e) = -(2.00314..) |/i| 8/15 
it easy to check that the the sign of the variations 

r r /l ( 2 fa 



m\m a 



h 



+ 0{t') 



(16) 



of the mass ratios r a = Rem a /mi agrees with the McCoy- Wu scenario (Fig. 7). 
The above values of f a b c give [7] 



Im 777.4 — 



-173.747..) r 2 



(17) 



The value (|15jl refers to the normalization of e in which (e(a;)e(0)) — > \x\ 2 as 
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Im m\ ~ (-49.8217..) r 2 , (18) 
and in turn the decay widths and lifetimes 

„ Imm 2 . 1 , , 

Ta = a ~ , t a = — 19 

m a T a 
of the particles and A$. The lifetime ratio 

lim — = 0.23326.. (20) 

T — t h 

is universal, as well as the branching ratio for A$, which decays at 47% into A\Ai and for the 
remaining fraction into A\A2- 

It appears from (|20|) that A^ lives more than four times longer than A4, somehow in contrast 
with the expectation inherited from accelerator physics that, in absence of symmetry obstruc- 
tions, heavier particles decay more rapidly. Notice, however, that in d dimensions the width for 
the decay A c — » A a Ab is 

r c ^ab^9 2 \fabc\ 2 ^d, (21) 

where g is the perturbative parameter, f abc the form factor and 

~ / 0~ 6 ( Pa +Pb~ Pe) (22) 

J Pa Pb m c 

the phase space (p = p a = ~Pb)- For fixed decay products, \p\ increases with m c and in d = 2 
suppresses the phase space, in contrast with what happens in d = 4. The results for the vertices 
f a b c indicate that in our case the dynamics further enhances the increase of t c with m c . 



4 Unstable particles in finite volume 

The issue of obtaining numerical checks of theoretical predictions for decay processes is made 
particularly interesting by the difficulty of characterizing unstable particles in the finite volume 
[22]. The problem is particularly relevant for lattice gauge theories. 

In two dimensions energy spectra can be obtained by numerical diagonalization of a truncated 
Hamiltonian on a cylinder geometry. The signature of particle decay on the cylinder is clear. 
At an integrable point, when the energy levels are plotted as a function of the circumference 
R of the cylinder, the line corresponding to a particle above threshold crosses infinitely many 
levels which belong to the continuum when R = 00 (Fig. 8). Once integrability is broken, this 
line "disappears" through a removal of level crossings and a reshaping of the lines associated to 
stable excitations (Fig. 9). 

One way of extracting the decay width from such a spectrum is the following [22]. Consider 
energy levels corresponding to states with two particles of mass m and momenta p and —p, 
sufficiently close to the threshold E = 2m that the particles can only scatter elastically, with 
scattering amplitude S(p) = exp iS(p). On the circle, p is quantized by the condition e lRp S(p) = 
1, or equivalently 

Rp + §(p) = 2vrn , (23) 
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Figure 8: First eight energy levels of the finite volume Hamiltonian of Ising field theory at 
r = as functions of r = m\R (from Ref. [4]). At r = 40, starting from the bottom, the levels 
are identified as the ground state, the first three particle states A\, A2 and ^3, three scattering 
states A±A±, the particle above threshold A4. Crossings between the line associated to the latter 
and the scattering states are visible around r = 18, r = 25 and r = 36. 
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Figure 9: First eight energy levels of the finite volume Hamiltonian of Ising field theory slightly 
away from r = (from Ref. [5]). Observe the splitting of the crossings pointed out in the 
previous figure. 
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I fcab | 


Exact 


Numerical 


1/411 1 
1/511 1 


36.730.. 
19.163.. 


36.5(3) 
19.5(9) 



Table 3: The values of the three-particle vertices f a b c 
compared with the exact predictions of [7J. 



fabc/ (e) obtained numerically in 



with n labelling the states. Hence 5{p) can be determined from the measure of 



E{R) = 2vVT 



(24) 



A narrow resonance of mass m c and width T can then be fitted through the Breit-Wigner formula 

E — m c — iT/2 



6(p) = 5o(p) + S B w(p) , 



$bw = — i In 



(25) 



E-m c + iT/2 ' 

5q being a smooth background. 

An alternative method was proposed in [7J. At an integrable point consider a particle with 
mass m c > m a + m&. The states 



II) 



\Mp)M-p)) = |2> 



(26) 



are degenerate on the cylinder at a crossing point R*. When integrability is broken by a per- 
turbation 

V = X / ds*(z), (27) 
Jo 

the energy splitting at R* is, to lowest order in perturbaton theory, 



AE = y/(Vu -V 22 ) 2 + MV 12 \ 2 , 



(28) 



where Vij = (i\V\j). Choosing a crossing point at R* large enough, the V^j's are well approxi- 
mated by the infinite volume matrix elements. In this way the decay vertices V12 can be obtained 
from measures of the energy splittings. 

Both methods have been used in [23] to measure the decay widths of the first two particles 
above threshold in Ising field theory close to the critical temperature from numerical diagonal- 
ization of a truncated Hamiltonian, directly in the continuum limit [M] . The results are shown 
in Table 3 together with the predictions of form factor perturbation theory discussed in the pre- 
vious section. It seems obvious that any numerical method designed to measure decay widths 
in lattice models should first of all be able to recover these predictions for the two-dimensional 
Ising model. 
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